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Recent hydrodynamical calculations for Au+Au collisions at ^snn = 130 GeV en- 
ergy are reviewed, and the initial conditions of hydrodynamical evolution necessary to 
reproduce experimental data are discussed. 

1. Introduction 

One of the first measurements at RHIC was the elliptic anisotropy of the particles pro- 
duced QTJ| . It was found to be as large as the hydrodynamical prediction |2| and subsequent 
hydrodynamical calculations achieved good fits to its centrality and pt dependence [|,(|,[5| . 
In the following I review how well hydrodynamical models can fit the data published after 
those early papers and what constraints the data has set on the initial state of hydrody- 
namical evolution. 

2. Comparison with the data 

As an example of hydrodynamical results, I mostly quote the work of Kolb and Heinz || . 
This is not meant to imply that other calculations are of lesser importance, but rather to 
emphasize how much of the data can be explained using one set of parameters. 

2.1. Transverse momentum spectra 

The conventional method to initialize hydrodynamical calculations is to fix the initial 
densities (energy or entropy and net baryon density) to reproduce the observed particle 
abundancies and to choose freeze-out temperature to reproduce the observed slopes of 
the spectra. This approach was also applied in || where the initial state was fixed to 
reproduce the charged particle multiplicity as a function of centrality and the freeze-out 
temperature was chosen to fit the central collision data. Therefore the good fit to pion 
and antiproton spectra in the most central collisions shown in the upper left panel of 
fig. is not surprising. On the other hand, the slope of the spectra as a function of 
centrality are predictions and, as can be seen, the agreement with data is impressive. For 
antiprotons the prediction lies within the experimental errors up to px = 3 GeV/c. Only 
for very peripheral collisions, b > 10 fm, do the data show a significant excess of particles 
at pt > 1.5 GeV/c. In hydrodynamical description the observed excess of antiprotons 
over pions at pr^,2 GeV/c is simply a consequence of strong transverse flow. 

In this calculation local chemical equilibrium is assumed to hold until kinetic freeze-out. 
Since the favoured kinetic freeze-out temperature, Tf = 130 MeV, is much smaller that 
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Figure 1. Charged pion, antiproton and positive kaon spectra from central (upper left 
panel) and semi-central to peripheral (other panels) collisions from H. The data were 
taken by the PHENIX § and STAR g] collaborations. 



the chemical freeze-out temperature, T c h = 165-175 MeV given by thermal models 0, it is 
not possible to obtain correct proton and antiproton yields simultaneously. In |J this was 
circumvented by using pion yields at Tf = 130 MeV as they were, but scaling by hand all 
the other particle yields to their calculated values at hadronization temperature T c = 165 
MeV. The results obtained in this way are very similar to those of || where the hadronic 
phase is described by a transport model and separate chemical and kinetic freeze-outs are 
included. Thus one can consider the results of || to be reasonable approximations. 

2.2. Elliptic anisotropy 

Since the initial particle production is azimuthally symmetric, azimuthal anisotropy 
of the final particle distributions is a signal of rescatterings among the particles pro- 
duced. More frequent rescattering can be expected to lead to a larger anisotropy, and 
since hydrodynamics assumes zero mean free path and thus an infinite scattering rate, 
it provides an upper limit to observable anisotropies. Anisotropy is quantified by mea- 
suring the harmonic coefficients v n (y,pT] b) of a Fourier expansion in <ft p of the measured 
hadron spectrum dN/ (dypT dpr d4> p ) [K|. Anisotropy characterized by a non-zero second 
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Figure 2. The elliptic anisotropy coefficient v 2 (px) of all charged particles (left) and 
identified pions and protons (right) at midrapidity in minimum bias collisions from M. 
The data are from the STAR |H, |T3| , |T4|1 and PHENIX |15|] collaborations. The curves 
correspond to calculations using an equation of state with (Q) or without (H) a phase 
transition and (in the right panel) three different freeze-out temperatures (Tf = 128 MeV 
(dash-dotted), 130 MeV (solid) and 134 MeV (dashed)). 



coefficient, v 2 , is called elliptic anisotropy or elliptic flow [] I 



The coefficient v 2 as a function of transverse momentum px in minimum bias collisions is 
shown in fig. [| The calculation was done using an equation of state with a first order phase 
transition (EOS Q) and without a phase transition (EOS H). The observed anisotropy of 
charged hadrons shown in the left panel is seen to reach the hydrodynamical values up to 
Pt < 2 GeV. Above px = 2 GeV the hydrodynamical result keeps increasing monotonically 
whereas the data saturates, indicating incomplete thermalization of the high momentum 
particles. The pion and proton anisotropies depicted in the right panel also show the 
hydrodynamically predicted mass dependence: at low values of transverse momentum, 
the heavier the particle, the smaller its v 2 [§J. Interestingly at this conference the pion 
anisotropy was shown to deviate from hydrodynamical predictions around px = 1.5-2 
GeV, whereas the proton anisotropy reaches hydrodynamical values even at px ~ 2.5 
GeV [|T2[]. As can be seen the anisotropy of hadrons or pions is quite insensitive to the 
phase transition, but the anisotropy of protons shows a clear dependence: the curve 
calculated assuming a phase transition is very close to the data, but an EoS without a 
phase transition leads to too large a proton anisotropy at low px- 

The fit to the proton data shown in fig. |2] is worse than achieved in earlier calculations @, 
nj. This is due to a different freeze-out temperature, which was Tf = 120 MeV in j|, but 



was chosen to be Tf = 130 MeV in || to achieve a better fit to the px spectra. On the 
other hand, the hydro + cascade approach of reproduces well both the px spectra and 
differential anisotropy v 2 (px) of pions and protons without such ambiguity. Whether this 
points to a weakness in the ordinary Cooper-Frye description of freeze-out or is a result 
of differences in the initialization of the calculations remains to be seen. 

In fig. |] the centrality dependence of the anisotropy coefficient v 2 is shown. The cal- 




Figure 3. Centrality dependence of elliptic 
anisotropy of charged particles at midra- 
pidity. The data are analysed using 4th 
order cumulant method by the STAR col- 



laboration [|16 



Figure 4. Pseudorapidity dependence of el- 
liptic anisotropy of charged particles in min- 
imum bias collisions from ||18|| . The prelim- 
inary data are from [|19j. The curves cor- 



respond to two different parametrizations of 
the initial state. 



culation is similar to that in || and the data is analysed by the STAR collaboration [16 



using the 4th order cumulant method [fT? |. The fit to the data is slightly worse than 
shown before |§ since the 4th order method leads to smaller values of v<i than previously 
published [[[]. The new values are, however, within the systematic errors of the older 
analysis. The present data is below the hydrodynamical limit in all but the most central 
collisions, but the discrepancy is significant only for large impact parameters, b > 6 fm. 

In the calculations mentioned so far the expansion is assumed to be boost invariant. If 
the assumption of boost invariance is relaxed and expansion in all three dimensions is done 
numerically, the centrality and px dependence of the anisotropy at midrapidity is similar 
to that in the boost invariant case |Ti| . When compared to the excellent fit to data shown 



in fig. 0, the pseudorapidity dependence of elliptic anisotropy shown in fig. f| |18[ looks 
less satisfactory. The data reaches the hydrodynamical value only around midrapidity. 
On the other hand, even this result reproduces the data within a window of one to two 
units of pseudorapidity. This region already contains most of the particles produced. It is 
worth remembering that anisotropy in hydrodynamical models depends strongly on the 



initial shape of the system. The initialization used in is relatively simple and a more 
sophisticated initialization may lead to better fit to the data. Therefore it is premature 
to conclude, based on this data and calculation alone, that thermalization is reached only 
at midrapidity. 

2.3. Two-particle interferometry 

The analysis of two-particle momentum correlations known as HBT interferometry 
provides a method to study the space-time structure of the emitting source |p0| . It 
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Figure 5. The HBT radii -R S id e and R out from [|23| . The data are from [25|. The curves 
correspond to equations of state with a first order phase transition but different latent 
heat (CIII and CII) or with a smooth crossover between the phases (CI). 



has been predicted that a first order phase transition would lead to unusually large HBT 
radii [21]]. However, comparisons of calculations |6,22,23,j24| with data p5|j26[| have lead to 
the so-called HBT-puzzle: All calculations give a ratio of HBT radii i? out /-R S ide larger than 
one, but the experimental value is of order one. As shown in fig. [5] |23|, hydrodynamics 
usually leads to a too small sideward radius -R S id e and to a too large outward radius -R ou t- 
This is usually interpreted to mean that the system expands more and its lifetime is 
shorter than given by hydrodynamics. 

So far it looks like it is possible to fit either R S id e or R out , but not both simultane- 
ously. A good fit to -R S ide was achieved in |22] where the hadronic phase was described 
using a hadronic cascade. This naturally leads to a spatial freeze-out distribution with 
a larger average size in the sideward direction and thus to larger -R S id e - An acceptable 
fit was achieved also in p3 where the initial size of the system was larger than in other 
calculations, but the outward radius -R out is too large in both of these works. In || it 
was shown that if freeze-out takes place immediately after hadronization, Tj = T c = 165 
MeV, -R out is reproduced, but -R S id e is even smaller than when the freeze-out temperature 
is Tf = 130 MeV. The effect of such a high freeze-out temparture on single particle spectra 
or anisotropics was not checked either. 

One interesting detail shown in fig. |5] is that the radii are closest to the data when an 
equation of state with a smooth crossover is used. One can thus claim that HBT radii 
favour an equation of state with crossover but, as argued above, proton v 2 results favour 
an equation of state with a first order phase transition. The explanation to this seemingly 
contradictory behaviour, as well as to the entire HBT-puzzle, is still unknown at present. 
For the most recent developments see p7 |. 



6 




2 4 6 8 10 
r (fm) 



Figure 6. Initial energy density distributions in the transverse plane scaled to r = 1 fm/c 
from UJ (Morita et al.) and Q (Eskola et al). 



3. Initial conditions 

It is well known that even if the hydrodynamical evolution is constrained to reproduce 
the observed single particle spectra, there still is considerable freedom in choosing the 
initial state. One degree of freedom which complicates the comparison of the initial state 
of different hydrodynamical calculations is the initial time, tq, of the calculation. Its effect 



is easily seen in the Bjorken estimate [28,29] for the initial energy density, 



dE T 1 

where a smaller initial time leads to a larger initial energy density and vice versa. Thus 
a meaningful comparison of initial states from different calculations requires scaling the 
initial densities to the same initial time r. This is usually done by assuming one dimen- 
sional Bjorken expansion with an ideal gas equation of state from the initial time tq of a 
particular calculation to a chosen common time r: e = eo( r o/ r ) 4 ^ 3 - 

However, even if the difference in initial time is taken into account, the final single 
particle spectra allow very different initial energy distributions. As an example of this, 
initial energy distributions used in recent calculations by Morita et al. [^4] and Eskola et 



al. 1 30] are shown in fig. |(| The former is a simple generalization of the one dimensional 



Bjorken model to transversely expanding central collisions: a flat energy density distri- 
bution with gaussian smearing near the edges. The latter is based on a pQCD saturation 



calculation [ET| of the transverse energy of minijets produced in the primary collisions 
which leads to a very peaked distribution of energy density. In this approach the net 
baryon density is also obtained from pQCD calculation and — after fixing the equation 
of state and equating the initial time with the formation time, r = l/p sa ,t — the only 
free parameter is the freeze-out temperature |32|. 



7 



c 

O 

i 

15 



5 



1000 




100 




10 


#\ 


1 




0.1 


•v 

A. V 


0.01 


0.001 




0.0001 




1e-005 





ir 

PHENIXtt 

r(xo.6) 

PHENIX/T 
p(x3.5) 
PHENIXp 



0.5 1 1.5 2 2.5 3 3.5 4 

m r [GeV/ c 2] 




Pt [GeV] 

Figure 7. The pr distributions of negative hadrons calculated by Morita et al. 



(left) 



and Eskola et al. [pOj (right) compared to data by PHENIX collaboration [p,33 



Even though the initial states in these two calculations are very different, the final single 
particle distributions shown in fig. [7] are surprisingly similar and reproduce the slopes 
equally well! The particle yields must again be treated with caution. In both of these 
works chemical equilibrium is assumed to hold until kinetic freeze-out. In [|24|] the kinetic 
freeze-out temperature Tf = 125 MeV is well below the chemical freeze-out temperatures 
given by thermal models. Thus the authors note that their calculation reproduces the 
observed slopes of the spectra, but to reach the yields, they must scale them by hand. On 
the other hand, the authors of |K| employ a considerably high freeze-out temperature, 
Tf = 150 MeV, which allows them to reproduce also the particle abundancies. 

Unfortunately other observables are not helpful in differentiating these two approaches 
either. Both of the calculations were done for central collisions where the elliptic anisotropy 
is zero. The HBT radii are sensitive to the build-up of flow and the lifetime of the system 
and therefore to the initial distributions, but unfortunately HBT radii were not calculated 
in |Kj . As an estimate one can use the HBT radii for a high freeze-out temperature shown 
in || instead. The radii in || and [24| are indeed different, but neither of them repro- 
duces the data: the first leads to a correct R out but too small a -R S idc, whereas the latter 
reproduces -R S id e acceptably but R out is too large. Another observable which is sensitive 
to the maximum temperature reached in heavy ion collisions is the spectrum of direct 
photons, but there is no photon data available yet, nor was the yield calculated in these 
papers. 

Even if the initial energy density distribution is only weakly constrained by the data, 
the initial energy per unit rapidity is almost fixed. In Table |l| the initial states of recent 
hydrodynamical calculations are characterized by scaling to r = 1 fm/c and calculating 
the average energy density in the transverse plane. The peak values of the initial energy 
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Table 1 

Average initial energy density in the transverse plane scaled to a time r = 1 fm/c, initial 
time, and freeze-out temperature in recent hydrodynamical calculations. 





(e)(l fm/c) 


t [fm/c 


Tf [MeV] 


Kolb & Heinz g 


| 5.4 GeV/fm 3 


0.6 


130 


Hirano |18| 


5.8 GeV/fm 3 


0.6 


137 


Teaney et al. || 


~ 6 GeV/fm 3 


1.0 


N/A 


Morita et al. [24 


| 3.9 GeV/fm 3 


1.0 


125 


Eskola et al. 30j 


6.5 GeV/fm 3 


0.19 


150-160 



density have a huge spread from 6 to 160 GeV/fm 3 in these works, but the average values 
at 1 fm/c have a much smaller spread from 3.9 to 6.5 GeV/fm 3 . This spread is partially 
due to the different values of the effective transverse area, irR 2 . If this is taken into 
account, the average energy density in all works is 5-6 GeV/fm 3 , very similar to the the 
Bjorken estimate 4.6 GeV/fm 3 f29[. 

The initial times of hydrodynamic evolution employed in these calculations have a 
relatively large spread from 0.2 to 1 fm/c. However, in all calculations the initial time 
is shorter than the perturbatively estimated thermalization time of a few fm/c ||34j| . The 
basic argument in favour of a short initial time is based on elliptic anisotropy; the estimate 
by Kolb et al. [Q] showed that if there is 1 or 2 fm/c delay in thermalization, the value 
of t>2 decreases by 10% or 25 %, respectively, which leads to values of v 2 which are below 
the data. 

There is no similar consensus about the freeze-out temperature at RHIC. The single 
particle spectra in central collisions can be fitted equally well (fig. 0) if freeze-out takes 
place almost immediately after hadronization at Tf = 150-160 MeV []30| or at Tf = 125 
MeV [p4j| , roughly at the same temperature as at the SPS. The elliptic anisotropy of 



protons provides an additional constraint on the freeze-out temperature. So far it seems 
to favour a relatively low freeze-out temperature of Tf < 130 MeV, but whether it is 
possible to fit the elliptic anisotropy data using same kind of initial state and freeze-out 
temperature as in [|30| remains to be seen. While discussing freeze-out temperature at 
RHIC it is also worth remembering that in the hybrid model of [||], where the plasma 
phase is described using hydrodynamics and the hadron phase using hadronic cascade 
(RQMD), freeze-out temperature is not well defined. 



4. Summary 

Hydrodynamical models have been very successful in explaining the single particle 
RHIC data at low pr- The pt spectra and anisotropies in central and semicentral collisions 
are well reproduced for px < 1.5-2 GeV, and the p/n ratio at pt ~ 2 GeV/c has a simple 
explanation due to flow. Especially impressive has been how hydrodynamics is able to 
create simultaneously an elliptic anisotropy of negative hadrons which is large enough and 
an anisotropy of protons which is small enough to fit the data. If one considers solely this 
data the collision system behaves like a thermal system. 
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However, the reproduction of the HBT radii has been unsuccessful so far. It is unclear 
whether one should refine the final freeze-out process, the hadronization process, or the 
initial state to reach an acceptable description of the data. Especially puzzling is the 
fact that the HBT radii seem to favour a relatively stiff equation of state with a crossover 
phase transition, whereas the elliptic anisotropy of protons favours a soft equation of state 
with a first order phase transition. 

The details of the initial state required to fit the data are not yet completely fixed. The 
common features of all calculations are that the collision system thermalizes rapidly - 
initial time is tq < 1 fm/c — and that the average initial energy density is well above the 
critical energy density. 
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